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Abstract 

A recent Monte Carlo study of quenched QCD showed that the chiral conden- 
sate is non-vanishing above T c in the phase where the average of the Polyakov 
loop P is complex. We show how this is related to the dependence of the spec- 
trum of the Dirac operator on the boundary conditions in Euclidean time. We 
use a random matrix model to calculate the density of small eigenvalues and 
the chiral condensate as a function of arg P. The chiral symmetry is restored 
in the argP = 2tt/2> phase at a higher T than in the argP = phase. In 
the phase arg P = tt of the SU (2) gauge theory the chiral condensate stays 
nonzero for all T. 
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I. INTRODUCTION 



The relation between chiral symmetry breaking and confinement in gauge theories has 
been a subject of continuous interest. The picture of chiral symmetry breaking suggested by 
Banks and Casher relates the chiral condensate to the average density of small eigenvalues 
of the Dirac operator in the fluctuating background of the gauge field. The appearance of 
small eigenvalues in turn can be related to the fact that the QCD vacuum is populated by 
instanton configurations [[|. 

The pure QCD without quarks undergoes a phase transition at some temperature T c . It 
is interesting to see whether this transition affects the behavior of small eigenvalues of the 
Dirac operator in such a way that the chiral symmetry is restored. This motivates a Monte 
Carlo study of QCD with quenched quarks. In this paper we analyze the recent result || 
of such a study. Unlike previous similar works ||] the authors looked at the value of the 
chiral condensate in each of the Z3 symmetric directions in the complex plane in which the 
Polyakov loop acquires a nonzero expectation value above T c . It was found that: 

(a) The chiral condensate (if^ip) vanishes above T c in the phase with argP = 0. 

(b) In the phase with the argP = 2ti/3 the (ipip) is nonzero for some range of T > T c ! 

The last point (b) contradicts a naive expectation that the vanishing of the (jpip) has 
something to do with the rearrangement of the QCD equilibrium state at T c , e.g., a sup- 
pression of the instanton configurations in the deconfined phase. In this paper we show that 
this result can indeed be explained very naturally if one takes into account only the effect 
of the boundary conditions on the quark fields in the Euclidean time. Nevertheless the fact 
(a) that the vanishing of (ipip) coincides with T c tells us that the rearrangement of the QCD 
equilibrium state at T c also plays a role. 

II. POLYAKOV LOOP AND THE SPECTRUM OF THE DIRAC OPERATOR 

The partition function of the QCD without quarks can be written as a Euclidean path 
integral: 

Z = J VA^exp{~S[A]}; S = J d 3 xF^F, u . (1) 

An order parameter distinguishing the confined and the deconfined phases is the expectation 
value of the Polyakov loop: 
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P(x) = i Tr Pexp lw dtA Q (t,x) 



(2) 



The value of P is invariant under a gauge transformation periodic in the Euclidean time: 
U(l/T) = U(0). However, the action S is also invariant under aperiodic gauge transforma- 
tions: 



U{l/T)=zU(0), 



(3) 



with z = Jexp{i27r/3}, which do not upset the boundary conditions on because z com- 
mutes with all matrices in the gauge group. The value of the Polyakov loop acquires a phase 
exp{i27r/3} after such a transformation. This corresponds to a Z 3 symmetry transformation 
in an effective theory of the Polyakov loops. 

At small temperatures the expectation value of the Polyakov loop is zero and the Z3 
symmetry is exact. Above T c the Polyakov loop acquires a non- vanishing expectation value 
which can be related to the liberation of color [HI. 
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(a) (b) 

FIG. 1. A schematic contour plot of the distribution of the averaged value of the Polyakov 
loop: (a) — below T c , (b) — above T c . Only 2 isolines are drawn: one at the level of the saddle 
point between the peaks and another in between the heights of the peaks. 

To analyze the results of the Monte Carlo simulation it is helpful to consider a probability 
distribution of the Polyakov loop in a large but finite volume. It can be calculated as a 
histogram ||: For each configuration one finds the averaged (over x) P. The configuration 
then contributes a unit to the height of the histogram at this P. The qualitative behavior 
near T c of this distribution is shown in Fig. |]. There are 4 peaks. At low T the peak at 
P = is taller. Most of the configurations have very small P near the origin. At T = T c the 
3 degenerate peaks at nonzero P are (roughly) of the same height as the peak at the origin. 
At higher T the Z 3 peaks are taller. 
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The chiral condensate is measured by averaging the inverse of the Dirac operator over 
the gauge configurations. At high T one has a choice of averaging over ensembles of config- 
urations clustered near each of the Z 3 peaks. We refer to these ensembles as Z 3 phases. It 
was found in || that the result is different in the phases argP = and argP = 27r/3. 

Due to the Z3 symmetry an ensemble of gauge configurations in the phase argP = 27r/3 
can be mapped onto the ensemble at argP = by an aperiodic gauge transformation 
The Z 3 phases of pure QCD are equivalent. However, the spectrum of the Dirac operator 
is not invariant under such a transformation.^ We conclude that different behavior of the 
chiral condensate in the two phases is entirely due to the effect of the aperiodic gauge 
transformation on the spectrum of the Dirac operator. 

The action of the gauge transformation on the fermion fields changes the boundary 
conditions in the Euclidean time direction from antiperiodic: 

V(l/T,aj) = -^(0 >a j) (4) 

to twisted: 

iP(l/T,x) = -zip(0,x). (5) 

The eigenvalue problem for the Dirac operator is defined by the differential equation: 

(ip+m)^ = Xip, (6) 

with functions ip{x) satisfying given b.c. The b.c. in this problem are different in the phases 
with argP = and argP = 27r/3. 

Moreover, qualitatively one can see that the effect of a nontrivial Z 3 phase on the b.c. 
tends to lower the eigenvalues, because of less twisting required of ip- It is transparent if 
you consider as an example the free Dirac operator, i.e. = 0. In this case, the spectrum 
in the phase argP = is given by: A 2 = k 2 + ((2n + l)7rT) 2 . While in the argP = 27r/3 
phase it becomes: A 2 = k 2 + ((2n + l/3)7rT) 2 . The smallest eigenvalue moves from A = nT 
to A = ttT/3. 



1 A simple and well known manifestation of this fact is the lifting of the degeneracy between the 
Z% minima of the effective potential by the fermion determinant contribution. 
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III. RANDOM MATRIX MODEL 



To find the dependence of the spectrum on argP for arbitrary is a very difficult 
problem. Fortunately, the only information we need about the spectrum is the density of 
small eigenvalues p(A — * 0) which is related to the chiral condensate by Jl]: 

(W)=np{0). (7) 

To go further we need to make some approximations. As our first approximation we 
neglect the Euclidean time dependence of the gauge configurations responsible for small 
eigenvalues. This is not altogether unreasonable approximation at least at T of the order of 
T c and higher for the following two reasons. First, it is the infrared behavior of the gauge 
fields that is responsible for the small eigenvalues. It is known, that the IR behavior of QCD 
at high T is given by a 3-dimensional theory 0. Second, the correlation length at the first 
order phase transition in pure QCD may be large compared to 1/T — the extent in the 
time direction. For the SU(3) theory this is suggested by the fact that in the Potts model, 
to which this transition is related by universality arguments, the transition is very weakly 
first order. For the SU(2) theory the correlation length is infinite at T c and the theory is 
effectively 3-dimensional at T c . 

In this approximation we can separate t and x dependence in the eigenfunctions: 
ip(t, x) = J2 n e luJnt ip n {x). Thus (H) reduces to the eigenvalue problem for a t/iree-dimensional 
operator: 

[ij uj n + i-f A + jD + m}i) n (x) = \i[) n (x) (8) 

for each u n , where the gauge fields A , A depend only on x. The values of u n are determined 
by the b.c. in the Euclidean time: 

u n = ((2n + 1)tt - arg P)T. (9) 

To find the density of small eigenvalues (or, chiral condensate) for the operator (|8]) we 
use another approximation known as the (Gaussian) random matrix model ||. We choose 
a basis ip a of N functions^ (the same for all gauge field configurations and all uj n ) and write 
the t/iree-dimensional Dirac operator in this basis: 



2 In 4d the functions ip a can (but need not necessarily) be thought of as the exact zero modes of 
each of N individual instantons and X being the overlap matrix. In 3d a similar role can, perhaps, 
be played by periodic instantons or dyons. 
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(10) 

where the X is a complex N x N matrix. Then one can replace the averaging over the 
gauge fields by averaging over the elements of the matrix X. The main approximation of 
the model, which proves to be very good (and, probably, exact in some sense) in 4d ||, is 
that the probability distribution is Gaussian. In other words, we calculate: 

Z= [vXe W \-NC 2 TrXXnT\det N f ( ™ iX + iu} A (n) 
J 1 } n \ iXi + iu n m j 

This is actually an approximation to the full partition function of QCD with Nf flavors 
of dynamical fermions. We will be interested, however, in the limit Nf — > — quenched 
fermions. In fact, one can check that Nf cancels out in our final results anyway. Therefore, 
to simplify the formulae we shall put Nf = 1. It is also convenient to view each Matsubara 
mode as a flavor with a Matsubara mass u> n . From Z we can obtain the chiral condensate 
in the thermodynamic limit: 

T d 

(ibib) = lim lim — — \nZ. (12) 

m-+ov 3 ->ooN f V 3 dm K ' 

where V3 is the volume of the 3-space. 

The calculation [|]J§ of ([TT|) is similar in spirit to the mean field calculation in the 
large N sigma models (Hubbard-Stratonovich transformation). We introduce a new random 
matrix Y, the integral over which can be calculated in the saddle point approximation at 
large N. The first step is to rewrite the determinant as a Grassmann integral: 

Z = J VXV^>Vm exp [-NC 2 Tr XX f } exp * n ( %X + ^™ j * n | . (13) 

The spinors carry an implicit index a (a = 1, ... , N). Integration over X a ^ can be carried 
out which creates the four-fermion term: 

1 (14) 



NC 2 

This term can be rewritten using another auxiliary random matrix Y, 



Z = J VYV^V^ exp [-NC 2 Tr FF 1 "} exp | S *° ^ y ] ) ^ \ ' ( 15 ) 

where M is a diagonal matrix u n S n k- 
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The integral over Grassmann variables can be then performed: 

Z = J VY exp [-NC 2 Tr YY^] det J 

At large A" the saddle point equation for the matrix Y is: 

- C 2 Y + (Y + m) [(F f + m)(Y + m) + M(Y + m)- l M(Y + m) 



Y + m iM 
iM yt + m 



0. 



(16) 



(17) 



This has a solution in the form of a diagonal matrix with eigenvalues y n that satisfy (we put 
m=0 here): 

-y2„. Vn 



(18) 



This equation is analogous to the well known gap equation in the large N sigma models. A 
nontrivial solution y n = (l/C)^/l — uj 2 C 2 exists only for those n for which: 1 — uj 2 C 2 > 0. 
Using flTjD and ([1^) we find that the chiral condensate is given by: 



T I 

M) = lim lim — ( A^ Tr 

m— >0 Vs—kx V3 \ 



In the saddle point approximation at large A^ using (|18D we find: 



Y + m iM 
iM Y ] + m 



(19) 



(20) 



where n 3 = 2N/V3 and o; n = ((2n + l)7r — arg P)T . 

We see that the chiral condensate vanishes (continuously: (3 = 1/2 in this approximation 
when min n \u n \ > 1/C, in other words, when the temperature exceeds: 

1 



T x (argP) 



(21) 



C\n — argP| 

where argP is mod 2n.f\ This result is in accordance with our expectation that nonzero 
arg P favors chiral symmetry breaking. 

Our model at argP = is similar to a model suggested in where only the lowest 
Matsubara frequency was taken into account. This makes sense since we are interested in 
the behavior of small eigenvalues of the Dirac operator. Indeed, only the band ujq of the 
eigenvalues contributes in (|20|) when T is close to T x and < argP < 2tt. In our case 
keeping all Matsubara frequencies ensures periodicity in argP mod 2ir. 



3 In principle, C depends on T. For example, for T > T c one should expect C 1 ~ g 2 T. The 



equation (21) then defines T x self-consistently. It is not unreasonable to expect, though, that 
C ~ const near T c . 
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IV. DISCUSSION AND CONCLUSIONS 



We analyzed a recent Monte Carlo observation that the chiral condensate (ipip) does not 
vanish in some interval above T c in quenched QCD in the phase with argP = 2n/3. This 
fact suggests that the rearrangement of the QCD equilibrium state at T c (deconfinement) is 
not the only source of the restoration of the chiral symmetry in quenched QCD. 

We calculated (approximately) the chiral condensate at finite T taking into account the 
dependence of the spectrum of the Dirac operator on the b.c. for the quark fields in the 
Euclidean time. We found that the (ipip) vanishes continuously at a temperature T x (argP) 
which depends on the phase of the Polyakov loop (|2l[) . 

This suggests the following interpretation of the result 0. As we approach the decon- 
finement transition from T > T c side the chiral symmetry breaking in the argP = 2tx/2> 
phase occurs earlier, at a temperature T x (2ir/3) > T c . The chiral condensate is zero in the 
argP = phase because T x (0) < T x (27r/3). It turns out that T c > T x (0). Indeed, the 
result || suggests that T c is only about 7% lower than T x (27r/3), while in our approxima- 
tion T x (0) is about 1/3 of T x (27r/3). At T c the Z 3 phases become metastable and at lower 
temperatures the single P = phase dominates (see Fig. |l|). 

Being cavalier enough, one can conjecture that the effect of this transition on {ipip) can 
be taken into account very roughly by averaging over 3 values of arg P. This would lead to 
a relation between the condensate just below T c and the condensate in the argP = 27r/3 
phase just above T c : 

(W>)(T C -) « ^>(T c +,argP = 2vr/3), (22) 

which does not seem to be in a gross contradiction with || . In other words, 3 phases of the 
Polyakov loop "mix together" below T c . 

In a picture of the chiral phase transition in the quenched QCD suggested here the 
restoration of the chiral symmetry at finite T is mainly due to the effect of the b.c. in the 
Euclidean time suppressing small eigenvalues of the Dirac operator. The rearrangement of 
the QCD equilibrium state plays a secondary role. This goes in line with the fact that, 
as far as we know, the instanton density behaves rather smoothly across the phase tran- 
sition [Q,^n,|ll|] .R The effect of the b.c. in the argP = phase is stronger than in the 



4 We must emphasize here that the present discussion mainly concerns a model of QCD where 
the feedback of quarks onto the dynamics of gluons is absent. The idea that confinement is not a 
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argP = 27r/3 phase and in the low temperature P = phase. To put this into another per- 
spective one can note that antiperiodic boundary conditions for the quark fields are related 
to the Pauli exclusion principle.^ The fact that Pauli blocking leads to chiral symmetry 
restoration is rather well known in various models of particle and nuclear physics. 

Another interesting consequence of (|2lD is that in the SU(2) gauge theory: T x (n) = oo. 
This means that the (ipip) stays nonzero for all T in the arg P = n phase! 

This result can be argued for also beyond the approximation of the random matrix model. 
Recall that the theory at very high T > T c becomes effectively 3-dimensional. Moreover, 
fermions which would "decouple" in the argP = phase do not in the argP = n phase 
because the b.c. are periodic. Thus the value of (ipip) is determined by the 3-dimensional 
SU(2) gauge theory at zero temperature. The value of the (ipip) is nonzero in such a 
theory [Q. Moreover, it can be estimated to be of the order of g% purely on dimensional 
grounds. This means that the chiral condensate in the SU(2) theory in the argP = 7r phase 
at T ^> T c is of the order: 

(V^> ~ g 4 T 3 . (23) 
Recently there have been two attempts [0,0 to describe the data |3J using a Nambu- 



Jona-Lasinio model. With some tuning of parameters it is possible to achieve a reasonable 



fit to the data in the model ||15|| . The present paper emphasizes the role of the behavior 
of small eigenvalues as a function of argP and in this respect differs from a somewhat 
phenomenological approach of |I^| , |To| . This allows us to obtain new results, for example, fl2"T| ) 



or (gBJ). 

It would be interesting to test these results in a Monte Carlo simulation. For example, 
one can consider b.c. conditions on fermions with an arbitrary twist. This would add to our 
knowledge of the mechanism of the chiral symmetry restoration in QCD. 

Acknowledgments: The author would like to thank A. Kocic, J. Kogut, C. Korthals- 
Altes and R. Pisarski (who brought the result [[| to the author's attention) for useful dis- 
cussions and comments and the theory group of the Brookhaven National Laboratory for 
hospitality. The work was supported by the NSF grant PHY 92-00148. 



prerequisite of the chiral symmetry breaking has been also suggested by a study Q of the SU (2) 
theory with adjoint quarks. 

5 For a discussion of the relation between the Pauli principle and the decoupling of fermions from 
a dimensionally reduced theory at high T see |]13|| . 
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APPENDIX A: THE DENSITY OF EIGENVALUES 



The eigenvalue density p(A) for the Dirac operator can be found in a closed form in our 
approximation. Let us consider only the lowest band of eigenvalues with u = min n |a; n |. 
Contributions of the other bands are additive. To find p(A) we note that one can calculate 
(?/>^)(m) for arbitrary m. On the other hand, p(A) is the discontinuity of {ij)ij))(m) along 
the imaginary axis, in other words: 



Re {ipip)(i\ + e) = vrp(A). 



(Al) 



The {ijjijj)(m) is given by ( |TP| ) if we do not take m — > 0. The necessary eigenvalue y can 
be found from (|17D by solving a cubic equation. Finally, we obtain: 



7rp(A) = V3 
C 2 



A + g) 1 / 3 + (VA - q) 



,1/3 



(A2) 



where one should take the branch (— 1) 1//3 = —1, and 



5 = ^( 1+ ^_^), 



A = — 

27 



A 2 



;i _ tff -'-{I- 20cu 2 - 8£ 4 + 4A 2 cD 2 ) 



(A3) 



where A = AC and Cj = ujC. 

The shape of this distribution is shown in Fig. |2|. At = the p(A) has a semicircular 
shape. As lu increases the distribution develops a dip at A = and when uj > (1/C) it splits 
into two bands leaving p(A = 0) = 0. 




FIG. 2. The eigenvalue density of the Dirac operator for a single band uj as a function of cj. 
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